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Abstract 

We derive the off-shell scattering matrix for a spherical scatterer. The result 
obtained generalizes the off-on-shell matrix commonly used in the theory of 
scalar waves propagation in random media. 
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I. INTRODUCTION 



It was recently recognized that the internal structure of scatterers can significantly alter 
the picture of wave propagation and localization in random media |]|§. When the size of 
scatterers cannot be neglected, scattering is anisotropic and a "dwell time" for waves appears 
in addition to the time of flight between scattering |2|||. This substantially modifies such 
parameters describing wave propagation as the diffusion coefficient @H-p lj , transport mean 
free path [H^Ol and transport velocity |T2|-|l4||. The microscopic approach to the problem 



of the renormalization of the diffusion coefficient is based on the Bethe-Salpeter equation 
for the field-field correlation function. This equation and the solution obtained involve the 
transfer matrix for a single scatterer, t kk ,, where k, k' are outgoing and incident momenta, 
respectively, and uj is the frequency J[5| . This matrix contains all information relating to a 



scatterer's structure, its influence on the anisotropy of scattering and on internal resonances. 
A great deal of this information is ignored in the commonly used on-shell approximation for 
the transfer matrix where both momenta k and k' are taken on the "mass shell" , k 2 =uj 2 . The 
off-shell transfer matrix, originally arising in the Bethe-Salpeter equation, is not restricted 
by this condition. This matrix resolves the structure of a scattered field at any distance away 
from a scatterer, while the on-shell matrix gives the far- zone asymptote of a field only. The 
on-shell approximation formally enters the theory of wave transport in random media via the 
commonly accepted ^-function approximation for the imaginary part of the Green function, 
ImGfc ia , oc 5(u 2 — c 2 k 2 ). The 5 approximation confines all transfer matrices in the expression 
of the diffusion coefficient on the mass shell. However, the presence of derivatives of t^, 
does not allow one to set uj = k automatically. Even within the on-shell approximation, 
one has to first calculate the derivatives and then to allow k — > uj. Neglecting this fact can 
lead to unexpected results. For instance, numerical calculations of the influence of weak 
absorption within the on-shell approximation show that the diffusion coefficient becomes 



larger than its value in the medium without microstructure [16|. This effect disappears if 
the off-shell matrix is used |fTD|jTT[| . 

The off -on- shell matrix for a scatterer with an infinite permittivity was first used in Ref. 
]T|]. Such a matrix gives a far-zone asymptote of a scattered field when the source is located 
at a finite distance from a scatterer. This destroys the symmetry between the incident 
and scattered momenta which exists due to the reciprocal principal for the Green function. 
Analysis of the renormalization of the diffusion coefficient employing the off-on-shell transfer 
matrix of a scatterer with a finite permitivity fllQj|lll has shown a strong enhancement of the 
previously obtained |pJJ^| corrections to D and qualitatively agrees with experimental data 
[17| . Numerical agreement can be attained with the use of an exact off-shell matrix. In this 



paper we derive the transfer matrix of a dielectric sphere and analyze some of its general 
properties. 

II. TRANSFER MATRIX OF A DIELECTRIC SPHERE 

The transfer matrix of a single scatterer is originally introduced via the Green function 
of the wave equation 

{uj 2 [l+e(v)] + A}G r ^ = -5 r _ rl , (1) 
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where r, r' are location vectors of the point of observation and the source, respectively; e(r) 
has a constant value, e, inside the scatterer and vanishes outside, the wave speed in this 
equation and below is equal to unity. Eq. (1) must be solved under conditions that G rr , 
and its derivatives dG r r , / dr are continuous in the entire space outside the source |L8| . For a 
spherical scatterer, one can obtain G r r , by utilizing separation of the spherical coordinates 
and using proper boundary conditions at the surface of the scatterer 
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where a is the radius of the scatterer, R = |r— r'| is the distance between the point of 
observation and the source and Q 2 = u 2 (l + e). In these equations and below we omit the 
index / denoting P(x) = P l (x), J x = J l+ i/ 2 (. x ) an d H x = H^ 1 , 2 (x), where P t (x) are the Legendre 
polynomials, Ji + i/ 2 { x ) are t ne Bessel functions and Hp rl , 2 {x) are the Hankel functions of the 
first kind. The coefficients of these expansions, obtained by boundary matching, have the 
form: 
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Here we use a short notation of the following determinant based on any pair of cylindrical 
functions UJx) and V„(x) 
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By definition, the transfer matrix and the Green function are related by the equation, 
Ckk , = ^kk'^fc + ^kk'^fe-^fe') where Dfr = [k 2 — oo 2 ] is the wave propagator in empty space. 
Calculating the Fourier transform of the Green function, J r J r , (27r)~ 3 G r r , exp (— ikr+zk'r'), 
and using expansions of the spherical and plane waves over the spherical functions, Y lm , one 
can obtain 
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The integrals of the Bessel functions arising in Eq. (5) can be found in Ref. [[HJ 

a 



J ax J fix 



(6) 



The integrals containing Hankel functions diverge on the upper limit. However, since we deal 
with the retarded Green function, the frequency is shifted in the upper half-plane, u—nv-HO- 
This gives an exponentially decaying factor in H^ r when r — > oo and provides convergence of 
the integrals. Using recurrence relations between cylindrical functions one can show 



x J a xHp x dx — a 2_p 2 
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(7) 



Using Eqs. (3), (6), and (7) one can represent the off-shell transfer matrix of a dielectric 
sphere in the form 
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where L% = k 2 — uo 2 is the inverse wave propagator in empty space. This expression has 
an explicit symmetry, t kk , =t k , k , as it must for a spherical scatterer. Eq. (8) cannot be 
simplified any more in a general situation. 

The off-on-shell matrix can be obtained from Eq. (8) by taking the outgoing momentum 
to the mass shell, k = u;k, and using the facts that | Ju a J uc \ — and IJ^H^] = —2i/na 
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The transfer matrix given by Eq. (9) coincides with the result of Refs. 0,|TT|]. In the limit 
e — > oo (Q —>■ oo) it reproduces the result of Ref. [1]]. This transfer matrix defines the 
far-zone asymptote of the scattered field in the case that the point source is located at a 
finite distance from the scatterer. That destroys the symmetry between the incident and 
scattered momenta. 

For a source at an infinite distance we have to take the incident momentum, k', in Eq. 
(9) to the mass shell. This results in the well-known expression for the on-shell transfer 
matrix of a spherical dielectric scatterer |21J 



^kk' 



E 

1=0 



*(Z+l/2)P(kk')|J^ Q 

^\ H u,a J Qa\ 



(10) 
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III. RAYLEIGH LIMIT 



All expressions given by Eqs. (8-10) can be substantially simplified in the Rayleigh limit 
when the scatterer size is the smallest linear scale in a problem. Taking into account condi- 
tions um, Qa, ka, k'a <C 1 and using the asymptotes for cylindrical functions one can show 
that 



^2r(/ + 3/2)r(/ + 5/3)' 
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Substitution of these equations into Eq. (8) yields the transfer matrix for long waves 
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Due to the presence of two factorial functions in the denominator of Eq. (11), the main term 
in the sum is the term with 1 = 0. This shows that in the Rayleigh limit an isotropic S 
term dominates in the scattered field even at distances of the order of the wavelength from 
the scatterer. Calculating the isotropic term of the transfer matrix at the mass shell we 
reproduce the well-known Rayleigh scattering amplitude, t(u)— eoo 2 a 3 /67Y 2 . 

Eq. (11) contains three types of terms with different behavior at the mass shell. The 
first type of term remains finite when k, k' — > u, the second type vanishes as L£, and the 
third type vanishes as (L k ) 2 . The third type of term were usually missed in calculations, 
since the off-on-shell transfer matrix was used instead of the off-shell matrix. This leads to 
incorrect values of the partial derivatives dt kk /duj 2 \k=uj and dt kk /dk 2 \k =UJ in the expression 
for the diffusion coefficient. Using Eqs. (8), (11) enables one to obtain correct results for 
these derivatives. Below we obtain some general relations for dt kk /du 2 and dt kk /d\s. which 
allow one to avoid lengthy straightforward calculations. 



IV. DERIVATIVES OF THE T MATRIX 



'ID 



Calculating the derivative of the operator equation G = D+DtD with respect to u 2 and 
using the fact that dD / du 2 = D 2 one can obtain 



^- = D 2 + D 2 tD + DtD 2 + D^-D 



(12) 



On the other hand, differentiation of Eq. (1) yields 
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q^ = G (1+ S)G = (D+DtD) (l+S)(D+DtD) (13) 

Here the operators D and £ are diagonal in the momentum and coordinate representations, 
respectively, (k|D|k') = 5^ (k 2 — uo 2 ) 1 , (x|£|x') = e(x)5 xx /. From the perturbation series 
expansion of the scattering operator, t = oj 2 £ + u A £D£ + . . . , it follows that £Dt = tD£ = 
tjuj 2 —£. Using this relation to eliminate £ from Eq. (13) and comparing the result with 
Eq. (12), we obtain the operator equation 

m t+tfv+By. (Mi 



did 2 LJ 2 \ LU 



Using the explicit form for matrix elements of all operators in Eq. (14), we obtain an ex- 
pression for the derivative dt k]l , / du 2 

<^kk' _ ^kk' , [ 4-ui -iw q 2 dq 

''kq^qk' 9 / 9 9N 2 • \ L0 ) 



dcu 2 u 2 J kq qk V(g 2 -w 2 ) 5 

A similar formula for derivatives of with respect to wave vectors or their combina- 
tions can be obtained by differentiating the momentum representation for the perturbation 
expansion of the scattering operator, 

^kk' — ^ 2 ^kk' + ^ / £kqA 1 £qk'+ Cj6 / / ^kq-^q^qq'A^q'k' ■ ■ ■ ■ 
Jq JqJq' 



As a result we have 

^kk' 2<^kk' . 2 f ^kq . u> , , ,, , 

~dk +UJ J (q 2 -U 2 ) dk (16) 

Here £ kk / =£(k— k') is the Fourier transform of the dielectric function of a scatterer, e(x). 
In the case of a spherical scatterer with the constant permitivity £(k) has the form that is 
easy to differentiate £(q) = (sin qa—qa cos qa) ine/q 3 . 

Eqs. (15) and (16) are particular convenient when both momenta k and k' are taken on 
the mass shell. In this case on-shell and off-on-shell matrices only appear in the expressions 
for the derivatives and one does not have to differentiate the cumbersome Eq. (8) and then 
calculate the limit of k — > cu. 
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